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Abstract. Among the processes governing the energy bal- 
ance in the mesosphere and lower thermosphere (MLT), 
the quenching of CChfi-^) vibrational levels by collisions 
with O atoms plays an important role. However, there is 
a factor of 3-4 discrepancy between the laboratory mea- 
surements of the CO 2 -O quenching rate coefficient, kyj, 
and its value estimated from the atmospheric observa- 
tions. In this study, we retrieve kyj in the altitude region 
85-105 km from the coincident SABER/TIMED and Fort 
Collins sodium lidar observations by minimizing the dif- 
ference between measured and simulated broadband limb 
15 pm radiation. The averaged kyj value obtained in this 
work is 6.5+ 1.5 x 10 ~ 12 cm 3 s -1 that is close to other es- 
timates of this coefficient from the atmospheric observa- 
tions. However, the retrieved ky t also shows altitude de- 
pendence and varies from 5.5 + 1.1 x 10 “ 12 cm 3 s~* at 90km 
to 7.9+ 1.2x 10 - 12 cm 3 s _1 at 105km. Obtained results 
demonstrate the deficiency in current non-LTE modeling of 
the atmospheric 15 pm radiation, based on the application 
of the CO 2 -O quenching and excitation rates, which are 
linked by the detailed balance relation. We discuss the pos- 
sible model improvements, among them accounting for the 
interaction of the “non-thermal” oxygen atoms with CO 2 
molecules. 


1 Introduction 


Infrared emission in 15 pm CO 2 band (/15 |j m ) is the domi- 
nant cooling mechanism in the Earth’s mesosphere and lower 
thermosphere (MLT, see also Appendix A for the abbrevia- 
tions not explained in the text for readability’s sake) (Gordi- 
ets, 1976; Dickinson, 1984; Goody and Yung, 1989; Sharma 
and Wintersteiner, 1990). On Earth, the magnitude of the 
MLT cooling affects both the mesopause temperature and 
height; the stronger the cooling, the colder and higher is the 
mesopause (Bougher et ah, 1994). This process is also impor- 
tant for the energy budgets of Martian and, especially, Venu- 
sian atmospheres (Bougher et ah, 1999), where CO 2 cooling 
compensates for the EUV heating of the dayside upper atmo- 
sphere. The 1 15 jj m radiation is used to retrieve vertical tem- 
perature distributions T (z) in Earth’s atmosphere by a num- 
ber of satellite instruments: the CRISTA (Offermann et ah, 
1999), the SABER (Russell et al„ 1999), the MIPAS (Fis- 
cher et ah, 2008) and in the Martian atmosphere the MGS 
TES observations (Smith et ah, 2001; Feofilov et ah, 2012). 
It is generally accepted that the main mechanism linking the 
15 pm CO 2 atmospheric radiation to the heat reservoir (trans- 
lational degrees of freedom of atmospheric constituents) is 
the inelastic collision of CO 2 molecules with 0( 3 P) atom; 
first, atomic O excites the CO 2 bending vibrational mode 
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during the collision: 

C0 2 (v 2 ) + 0( 3 P) -> C0 2 (v 2 + 1) + 0( 3 P) (1) 

after which the excitation may be quenched either by an- 
other collision with some molecule or atom or by emis- 
sion of the radiation quantum: C0 2 ( u 2 + 1 ) — > C0 2 (t> 2 ) + 
hv { 667 cm -1 ), where i> 2 is the bending mode quantum num- 
ber. Both the cooling efficiency and /15 (j m strongly depend 
on the rate coefficient of the process (1) and on the atomic O 
volume mixing ratio (VMR). To be consistent with a gener- 
ally accepted way of describing this process, we will refer to 
the rate coefficient of the reaction inverse to Eq. (1) and will 
call it the “C0 2 -0 quenching rate coefficient” or Avt, where 
VT index stands for vibrational-translational type of interac- 
tion. Generally, it is assumed that the velocity distribution of 
atomic oxygen is Maxwellian and that the fine structure of 
atomic oxygen does not affect the process (1) and its inverse. 
First, we will use these assumptions that are typical for at- 
mospheric modeling and then will address their applicability 
in the discussion part of the work (Sect. 4). 

It is self-evident that both the calculation of radiative cool- 
ing/heating rates in C0 2 and the interpretation of measured 
1 15 (j m radiation require the best possible knowledge of the 
Ayr (see e.g., Garcia-Comas et al., 2008). However, despite 
the importance of ky t for the atmospheric applications, the 
values, obtained in the laboratory and retrieved by fitting 
the space observations, vary by a factor of 3-4 (see Table 1 
and Sect. 2 below for more details). In this work, we de- 
scribe the retrieving of the Avt from the coincidental space 
and lidar observations. For this purpose, we used an exten- 
sive data set provided by the SABER instrument aboard the 
TIMED satellite which contains, besides other information, 
vertical profiles of / 15 jj m (z) limb radiation, O(z), and, more 
recently, CO 2 (z) (Rezac, 201 1) VMRs. This dataset was sup- 
plemented with T(z) in 80-1 10 km altitude range measured 
by the Fort Collins lidar (40.6° N, 105.2° W). We show that 
the synergy of these two instruments enables one to retrieve 
kyj and study its behavior in the MLT. The comparison of 
the SABER temperature retrievals with the Fort Collins li- 
dar measurements has already been done by Remsberg et 
al. (2008) who found that these data sets agree consider- 
ing 50% error in quenching rate (Avt — (6.0 + 3.0) x 10 -12 
[cm 3 s -1 ]). In this work, we try to make a step forward and 
get into more details of the MLT physics by applying strin- 
gent overlapping criteria and by treating the individual al- 
titude layers of the overlapping region separately. The dif- 
ference of the present study from the work of Remsberg et 
al. (2008) is in the kyj retrieval algorithm: instead of compar- 
ing the temperatures retrieved with different kyj values, we 
perform the kyj(z) estimate by minimizing the differences 
between the measured and calculated /15 pim(z)- All the cal- 
culations presented in this work were carried out using the 
non-LTE ALI-ARMS code package (Kutepov et al., 1998; 
Gusev and Kutepov, 2003; Feofilov and Kutepov, 2012). The 
background for the non-LTE problem for the molecular gas 


and the review of kyj measurements and estimates is given 
in the next section. 


2 Non-LTE problem for the molecular gas in the 
atmosphere and Ay t rate coefficient 

Inelastic molecular collisions determine the population of 
molecular levels in the lower atmosphere. As a result, local 
thermodynamic equilibrium (LTE) exists where the popula- 
tions obey the Boltzmann law with the local kinetic temper- 
ature. In the MLT, the frequency of collisions is lower, and 
the vibrational level populations must be found taking into 
account all the processes, which populate and depopulate vi- 
brational levels: the absorption of atmospheric and solar ra- 
diation in the ro-vibrational bands, spontaneous and stim- 
ulated emission, chemical sources, vibrational-vibrational 
(V-V) and vibrational-translational (V-T) energy exchange 
processes. The altitude, above which the LTE approxima- 
tion is not applicable, depends on the relationship between 
these processes: in general, if the V-T processes dominate or 
the radiative processes are balanced (optically thick media), 
the vibrational levels populations are close to LTE. For the 
C0 2 (v 2 ) vibrational levels involved in formation of /15 jLi m , 
the non-LTE effects become significant above ~ 75-80 km 
altitude (Lopez-Puertas and Taylor, 2001; Kutepov et al., 
2006; Feofilov and Kutepov, 2012). 

The importance of kyj rate coefficient for the calcula- 
tion of C0 2 emission in the MLT was first discussed by 
Crutzen (1970). He suggested an estimate for this value with 
the upper limit of 3. Ox 10 -13 cm 3 s -1 . The first laboratory 
measurement of Avt was performed at high temperatures 
( T >2000 K) (Center, 1973) using shock tube technique. The 
extrapolations of these measurements to room temperatures 
by fitting the Landau-Teller expression (Taylor, 1974) pro- 
vided the values of about 2.4x 10 -14 cm 3 s -1 . The average 
value of this rate coefficient obtained in later studies has 
changed by two orders of magnitude, and since the middle 
of 1980-s, the kyj is accepted to be on the order of (1.0- 
10.0) x 10 -12 cm 3 s -1 (see Table 1). However, as we show 
below, these variations are still large both for the adequate 
estimation of the radiative budget of the MLT region and for 
temperature retrievals from 15 pm C0 2 radiation. There is 
a well-known discrepancy between the laboratory measure- 
ments of Avt and its retrieval from the atmospheric measure- 
ments. Generally, laboratory measurements provide low val- 
ues of Avt centered around 1.3 x 10 -12 cm 3 s -1 . Huestis et 
al. (2008), based on the analysis of experimental data and 
quantum-mechanical calculations, recommend using Avt = 
1.5 x 10 -12 cm 3 s -1 . However, applying this value to the in- 
terpretation of /i5|j m measurements leads to overestimating 
the MLT temperatures. Instead, the values required for an ad- 
equate interpreting of atmospheric measurements performed 
in recent ~ 20 yr are usually about 5.5 x 10 -12 cm 3 s -1 with 
the exception of Avt = 1.5 x 10 -12 cm 3 s -1 retrieved by 
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Vollmann and Grossmann (1997) from the sounding rocket 
observations. 

To demonstrate the influence of kyj on the MLT area, 
we performed a sensitivity study for an average midlati- 
tude summer atmospheric scenario from MSIS-E-90 model 
(Hedin, 1991) using a standard set of V-V and V-T rate co- 
efficients (Shved et al., 1998; Kutepov et al., 2006) and the 
ky t- which was first set to 1.5 x 10 -12 cm 3 s -1 and then to 

6.0 x 10 _ 12 cm 3 s -1 . The results are presented in Fig. la-c. 
Figure la demonstrates the population of the first 1 ’ 2 -excitcd 
level shown as vibrational temperature (see the figure caption 
for the vibrational temperature definition). Figure lb shows 
the sensitivity of the total cooling/heating rate to the kyj 
change. After obtaining the non-LTE populations of all vibra- 
tional levels involved in the task, the broadband 1 15 jj m was 
simulated in line by line mode and the resulting spectrum was 
convolved with the “narrow” 15 pm SABER bandpass func- 
tion for two test runs. The corresponding change in the limb 
radiance is shown in Fig. lc. As one can see, the MLT area is 
sensitive to the kyj changes above ~ 85 km altitude, which 
is used as the lower limit for the ky t retrieval in this work. 
Below this level, the sensitivity rapidly decreases, and 85 km 
altitude can be considered as a lower limit for the kyj es- 
timates from /15 ^observations. The upper limit of 105 km 
for the kyj retrieval is defined by the fading of the signal 
strength with increasing altitude and by the uncertainties in 
the O (z) distributions. 


3 Retrieving ky x from the overlapping SABER and 
lidar measurements 

3.1 The ky t retrieval approach 

The general idea for the kyj rate coefficient retrieval from 
overlapping satellite and lidar measurements is in mini- 
mizing the difference between the measured and simulated 
It 5 pm by varying the kyj. The simulations are performed 
with the “reference” temperature profiles measured by the li- 
dar instrument, which are not affected by the uncertainties 
in the kyj coefficient. A similar approach was utilized by 
Feofilov et al. (2009) who used the H 2 O VMR profiles mea- 
sured by the ACE-FTS instrument (Bernath et al., 2005) as 
reference ones and estimated three rate coefficients neces- 
sary for the calculation of H 2 CXV 2 ) populations in the MLT. 
Retrieving the kyj from comparing the measured and sim- 
ulated /i 5 pm is somewhat more complicated because the 
CCLL^) populations depend not only on kyj but also on 
atomic oxygen concentration (or VMR), which contributes 
to uncertainties in retrieved kyj. The way of overcoming this 
limitation will be described below. First, let us consider a 
simplifed case of a single overlap for which everything ex- 
cluding kyj is known. Since calculations demonstrate mono- 
tonic dependence of CCL ( i-> 2 ) populations and limb radia- 
tion on kyj at all altitudes (see Fig. la, c in this work and 


Sect. 3.6.5. 1 in Lopez-Puertas and Taylor, 2001), the devia- 
tion f (kyj, z) = | ^meas (z) - fsimul (kyj , z ) | will have a single 
minimum at each altitude z and in the ideal case of noiseless 
signal the retrieved rate coefficient will be unique. Adding 
noise to the experimental radiation / meas (z) and adding un- 
certainties to calculated radiation hi mu i(kyj, z), which are 
linked with uncertainties in lidar temperatures and spatiotem- 
poral variability of the area, respectively, will blur the mini- 
mum of f {kyj, z) that will, finally, define the uncertainty for 
the kyj retrieval. Let us now consider the case, for which 
both atomic oxygen and kyj are not known. This exercise 
is necessary since even though the SABER retrieves individ- 
ual O(z) profiles (Smith et al., 2010), any offsets or errors 
in O(z) will propagate to kyj. However, this problem might 
be overcome if the average O av er(z) profile is known with 
a sufficient accuracy (from climatology, modeling or other 
measurements). In this case, one can search for a minimum 
of f (y, z) with respect to a new variable y — kyj x O over 
a large number of the SABER/lidar overlaps. At this stage, 
individual distributions O(z) are not needed. It is crucial to 
choose a grid for y in such a way that the following criteria 
are satisfied: (a) y variation range includes the y inm value that 
corresponds to an absolute minimum of £(z, y); (b) the grid 
step is fine enough to hit the minimum of £(z, y). When the 
minimum of f (z, y ) is found over a large number of over- 
laps, one can retrieve the optimal value of rate coefficient: 
kyj — /min/Oaver for each altitude point where O ave r is the 
average value of atomic oxygen VMR obtained either from 
the SABER or the other sources. At this point, the accuracy 
of Oaver(z) becomes important and will be discussed below. 
From the methodological point of view, if the dependence 
of radiance on y is strongly non-linear, and the y variabil- 
ity at each altitude is high, then the retrieval approach should 
be modified and the corresponding t;(kyj,z) should include 
the weights, which are proportional to the radiance responses 
at given y. However, as the analysis performed for this work 
shows, non-linear effects have little influence on the retrieved 
kyj for the conditions considered in this work. 

3.2 Using the Colorado State Sodium lidar temperature 
measurements for kyj retrievals 

In this study, we used /15 |j m (z), T(z), P(z), CC> 2 (z), and 
O(z) from the SABER VI. 07 database and coincidental T (z) 
measured by the Colorado State Sodium lidar described in 
details in (She et al., 2003). Briefly, the lidar is a two- 
beam system capable of simultaneous measurement of the 
mesopause region temperature and winds, day and night, 
weather permitting. This lidar was modified in 1999 in re- 
sponse to TIMED satellite objectives. The lidar setup can 
perform simultaneous measuring of the mesopause region 
Na density, temperature, zonal and meridional wind pro- 
files with both daytime and nighttime capability. The mea- 
surement precision of the lidar system for temperature and 
wind with 2 km spatial resolution and 1 h integration time 
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Table 1. Historical review of the /tvt!C 02 - 0 } quenching rate coefficient measurements and atmospheric retrievals at T = 300 K. 


fc VT {C 0 2 - 0 ) [crnV 1 ] 

Reference 

Comments 

3-30x10“ 14 

Crutzen (1970) 

First guess 

2.4x 10 “ 14 

Taylor (1974), Center (1973) 

Laboratory measurements 

5.0 x 10 “ 13 

Sharma and Nadille ( 1981 ) 

Atmospheric retrieval 

1.0 x 10“ 12 

Gordiets et al. (1982) 

Numerical experiment 

2 . 0 x 10 -13 

Kumer and James (1983) 

Atmospheric retrieval 

2.0 x 10“ 13 

Dickinson (1984); Allen (1980) 

Laboratory measurements 

5.2 x 10 “ 12 

Stair et al. (1985) 

Atmospheric retrieval 

3.5 x 10 “ 12 

Sharma (1987) 

Atmospheric retrieval 

3-9 x 10 “ 12 

Sharma and Wintersteiner (1990) 

Atmospheric retrieval 

1.5 x 10 “ 12 

Shved et al. (1991) 

Laboratory measurements 

3-6 x 10 “ 12 

Lopez-Puertas et al. (1992) 

Atmospheric retrieval 

1.3 x 10 “ 12 

Pollock et al. (1993) 

Laboratory measurements 

5.0x 10 “ 12 

Ratkowski et al. (1994) 

Atmospheric retrieval 

5.0 x 10 “ 13 

Lilenfeld (1994) 

Laboratory measurements 

1.5 x 10 “ 12 

Vollmann and Grossmann (1997) 

Atmospheric retrieval 

1.4 x 10 “ 12 

Khvorostovskaya et al. (2002) 

Laboratory measurements 

1.8 x 10“ 12 

Castle et al. (2006) 

Laboratory measurements 

6.0 x 10“ 12 

Gusev et al. (2006) 

Atmospheric retrieval 

1.5 x 10 “ 12 

Huestis et al. (2008) 

Recommended value 

1.3-2. 7 x 10 “ 12 

Castle et al. (2012) 

Laboratory measurements 



Fig. 1. The sensitivity of (a) CCHCOIO) main isotope vibrational level populations, (b) infrared cooling/heating rate in CCH bands, and 
(c) ^15 pm(~) to ^VT- F° r panels (a) and (b): lines with circles correspond to k\j j = 1.5 x 10 “ 12 cm 3 s -1 and lines with triangles corre- 
spond to fey T — 6 . Ox 10 “ 12 cm 3 s“ 3 . The non-LTE populations in panel (a) are presented as vibrational temperatures, which represent the 
vibrational level excitation against the ground level 0: n\/riQ = gi/go exp[— (E\ — Eo)/kT v ], where n oj, go,l, and £o,l are populations, 
degenerations, and energies of the ground state and first vibrational level, respectively. The atmospheric model used for this sensitivity study 
is the midlatitude summer profile of temperature, pressure, and atomic oxygen VMR from the MSIS-E-90 combined with the WACCM CO 2 
VMR profile for the corresponding conditions. The curves with squares/triangles in panel (c) show the /15 p m (z) responses for the cases 
when the kyj changes from 1.5 to 6.0 x 10 “ 12 cm 3 s -1 and varies within estimated uncertainties in Fig. 2c, correspondingly (see also the 
discussion in Sect. 4). 
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were estimated for each beam under nighttime fair sky con- 
ditions to be, respectively, 0.5 K and 1 .5 m s _ 1 at the Na peak 
(92km), and 5K and 15ms _1 at the edges (81 and 107km) 
of the Na layer. Depending on the purpose of the analysis, 
the temporal resolution may be made between 10 min and 
several hours. We have searched for the SABER/lidar simul- 
taneous common volume measurements in 2002-2005 using 
stringent criteria for time and space overlapping: A lat < 2°, 
Along < 2°, At < 10 min. Most of the profiles selected in 
way (85%) fall in 18-6 h local time interval. We substi- 
tuted the SABER T ( z ) in 80-1 10km altitude range with the 
corresponding lidar T(z) and hydrostatically adjusted P(z) 
to new T(z). Atomic O is changing in this period and us- 
ing individual O(z) is not reasonable; therefore, we used 
the y — ky txO variable discussed in Sect. 3.1. To reduce 
the number of runs with obviously incorrect y values, we 
used the following approach: at each altitude the grid on y 
was built using the available O(z) and 21 points for kyx in 
the (1.0-10) x 10 -12 cm 3 s -1 range with 5 x 10 -13 cm 3 s -1 
step. The correctness of the y grid selection was verified 
at each altitude and for each overlapping event by check- 
ing for the existence of the £(y, z) minimum. For each over- 
lapping event, the non-LTE populations of CO2 vibrational 
levels were found at all altitudes and /i5(j m was simulated 
in line by line mode and then convolved with the corre- 
sponding SABER bandpass function. This procedure was re- 
peated for all grid points of the y variable. Then at each al- 
titude point z we calculated the radiation difference ((y,z) 
(Fig. 2a). All f (y, z) curves up to 105 km (cyan curve) have 
a clear minimum £ m in()Az) that washes out for t,{y,z) de- 
pendencies above that altitude (not shown in Fig. 2a). This 
behavior is explained both by larger lidar T (z) and by larger 
/i5pm(z) uncertainties at higher altitudes. The values of y 
corresponding to £ m m(y,z ) form a separate scientific prod- 
uct y m in(z) that may be used in midlatitude atmospheric ap- 
plications for the /i5|j m calculations. The retrieved y lnm (z) 
values are 1.1 x 10“ 14 ; 5.7 x 10“ 14 ; 2.0x 10“ 13 ; 5.1 x 10“ 13 ; 
6.4 x 10 -13 [cm 3 s“‘ ] for z = 85; 90; 95; 100; 105 km, re- 
spectively. To obtain the kyx (z) profile (Fig. 2c) we di- 
vided the ]/ m j n (") profile by average atomic O VMR profile, 
Oaver(z), (Fig- 2b). The dashed lines in Fig. 2c represent stan- 
dard deviations estimated from input data uncertainties. 

4 Discussion 

Overall, the kyx(z) values shown in Fig. 2c fit well to the 
atmospheric retrievals: the averaged value of kyx is equal 
to 6.5 + 1.5 x 10 -12 cm 3 s -1 . However, Fig. 2c also demon- 
strates the altitudinal variability of kyx(z) that goes slightly 
beyond its uncertainties in 85-105 km altitude range. Obvi- 
ously, this variability does not imply that kyx rate coefficient 
depends on altitude. Fet us consider possible reasons for the 
observed kyx behavior. The retrieved kyx(z) depends on: 
(a) lidar T(z) in 80-1 10 km, (b) the SABER P(z) and T(z) 


below 80 km, (c) I [5 p m (z) (d) COzlz), (e) O(z), (f) CO2 non- 
LTE model. Offsets in any of these parameters will lead to 
offsets in the retrieved kyx(z). The quality of the SABER 
V 1 .07 product has been evaluated by Remsberg et al. (2008). 
As they show, in general, SABER VI. 07 temperatures are 1- 
3 K higher than lidar in the lower stratosphere and ~ 1-3 K 
lower than lidar in the upper stratosphere and lower meso- 
sphere. Assuming a 3 K positive temperature change in the 
lower stratosphere and 3 K negative temperature change in 
the upper stratosphere and lower mesosphere, we have es- 
timated / 1 5 jjm limb radiation change above 85 km. In this 
test performed for the mid-latitude atmospheric scenario, the 
pressure profile has been hydrostatically adjusted, the non- 
LTE populations have been found, and 1 15 (j m has been calcu- 
lated and compared to /15 |j m calculated for unperturbed pres- 
sure/temperature distribution. This test shows that changes in 
1 15 |j m above 85 km do not exceed 1 %. At the same time, the 
sensitivity of 1 15 jj m to kyx change within the error bars is 
much higher: 4% at 85 km, 8% at 90 km, 15% at 95 km, 
20% at 100 km, and 25 % at 105 km (see a solid curve with 
triangles in Fig. lc). Correspondingly, the propagation of the 
systematic error in pressure/temperature distributions below 
85 km to the area of our particular interest through hydrostat- 
ics is negligible. The main contributor to the 15 pm upwelling 
radiative flux at the mesospheric altitudes is the stratopause 
area where temperatures are high and the concentration of the 
emitters is still high. In this area, SABER temperatures are 
in reasonable agreement with other measurements (Rems- 
berg et al., 2008), and the uncertainty of the upwelling flux 
associated with temperature uncertainty can be estimated as 
+ 2 %. There are not too many CCMz) data sets for the MLT 
area. In Fig. 3a we show the average vertical distributions of 
CO2, which were calculated in WACCM model (Garcia et al., 
2007; Smith et al., 2011), retrieved from SABER V1.07 in 
(Rezac, 2011), and retrieved from the ACE-FTS occultation 
observations (Beagley et al., 2010). The distributions agree 
within standard deviation limits (shown only for ACE-FTS 
for readability of the plot, the others are of the same order of 
magnitude). We did not consider CRISTA CO2 distributions 
(Kaufmann et al., 2002) because there was no temporal over- 
lap with the data sets used in the present study: almost 10 yr 
have passed since the CRISTA-2 operation. 

The most critical component for the ky x(z) retrieval is 
Oaver(z). A detailed analysis of the SABER VI. 07 data 
shows that the O density is at least twice that from several 
other data sources (Smith et al., 2010 and references therein). 
However, reducing O ave r(z) by factor of two will mean in- 
creasing kyx(z) by the same factor, which will make it 8-10 
times larger than the laboratory measured values (line with 
squares in Fig. 3b). On the other hand, there is increasing ev- 
idence for the reliability of the O determined from SABER 
(Xu et al., 2012) and additional observations with similar 
high O VMRs (Sheese et al., 2011). The retrieval of O(z) 
in the SABER VI. 07 also depends on temperature. How- 
ever, as discussed by Smith et al. (2010), the contributions 
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Fig. 2. Estimating the optimal kyj from the overlapping SABER and lidar measurements: (a) deviations between calculated and measured 
/l 5 (j m at different altitudes with respect to a combined y value (see text). Each f (y, z) curve in this panel represents the average over 72 
individual deviations. Note the existence of f (y, z) minima at all heights up to 105 km (cyan curve); (b) solid line: average O(z) built for all 
overlapping events, dashed lines: standard deviations for O(z) estimated in Mlynczak et al. (2012); (c) solid line: A'vt(z) obtained as a result 
of dividing the individual minima found in the left panel by atomic O VMRs from the middle panel; dashed lines: standard deviations for 
F V t- 


of uncertainties in temperature to O(z) are small. Moreover, 
any appreciable errors in temperatures would upset the ex- 
cellent agreement between daytime and nighttime O (z) (see 
discussion in Smith et al., 2010). In this work, we have used 
the Of") uncertainties from (Mlynczak et al., 2012) for the 
kyj(z) error bars estimates shown in Fig. 2c. 

It is worth noting that the altitude dependence of the re- 
trieved kyj(z) follows the temperature distribution in the 
MLT. As the quantum-mechanical calculations show (de 
Lara-Castells et al., 2006), the vibrational quenching of 
CC> 2 (v 2 ) is significantly enhanced by spin-orbit interaction, 
the probability of which depends on energy, so positive tem- 
perature dependency should be expected. In the present study 
we already used a standard (7/300) ^ 2 temperature scaling 
of kyj, so one has to assume that the temperature depen- 
dency is stronger than 7 1 / 2 . However, the most recent lab- 
oratory measurements of ky t in the 1 42-490 K temperature 
range (Castle et al., 2012) do not confirm this prediction and 
demonstrate negative temperature dependency. In Fig. 3b we 
show kyj(z) estimates for this case. As one can see, kyj al- 
most reaches the value of 1 .2 x 10 “ 1 1 cm 3 s ~ 1 at 105 km. We 
believe that the contradiction between laboratory measure- 
ments and theoretical estimates of temperature dependency 
of ky t is just another hint for deficiency of the CO 2 non-LTE 
model with respect to CCL-O collisions, which will be dis- 
cussed further. 

The standard pumping term in the non-LTE model, which 
describes total production of CC> 2 (v 2 ) in the state with the 


number of bending mode quanta V 2 due to collisions with the 
Of 3 P) atoms has the form of 

Yv 2 = n O( 3 P){ W V 2 — 1^2— 1,V2 n\, 2 k v 2 ,y 2 — l} (2) 

where «o( 3 P) ' s the 0( 3 P) density, n vl -\ and n ,,, are the vi- 
brational states populations, and k V2 _\ V2 and k V2 jV2 _i are rate 
coefficients for one-quantum excitation and de-excitation, re- 
spectively. In current non-LTE models, including the one ap- 
plied in this study, it is usually assumed that k u ,_i v , = ko, 1 
and k V2 , V2 _i = A'l o- It follows from Huestis et al. (2008) that 
if the velocity distribution of 0( 3 P) atoms is Maxwellian and 
their fine structure is thermalized then the laboratory mea- 
sured ko. 1 and k\$ are linked by the detailed balance relation: 

^ 0,1 = £ 1,0 ■ — • e~ El/ kT ( 3 ) 

go 

where go and g \ are the statistical weights of the lower and 
upper vibrational states, respectively, E\ is the vibrational 
energy of the first V 2 vibrational level, k is the Boltzmann 
constant, and T is the local kinetic temperature. Sharma et 
al. (1994) showed that both above-mentioned conditions are 
valid for 0( 3 P) atoms in the Earth’s atmosphere up to at least 
400 km, which seems to justify usage of Eqs. (2) and (3) in 
the non-LTE models. However, as Balakrishnan et al. (1998) 
and Kharchenko et al. (2005) show, the non-thermal 0( 3 P) 
and ()( 1 D ) atoms are produced by O 2 and O 3 photolysis and 
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Fig. 3. (a) Average vertical distributions of CO 2 VMR distributions for 40.6 + 2° N in 2002-2009 for WACCM and SABER, and in 2008- 
2009 for ACE-FTS. Dashed blue lines denote standard deviation for the ACE-FTS CCD profile; (b) estimated vertical profile of Fvt(z) for 
O (z) divided by 2 and for a negative temperature dependency of &vt according to the latest laboratory measurements (Castle et al., 2012). 


Oj dissociative recombination reactions in the MLT. These 
“hot” atoms may serve as an additional source of CCMid) 
level excitation. Therefore, the expression (2) may need to 
be replaced by an expression like 


Yv 2 «0( 3 P) ' [(1 ttf) ■ |/Zy 2 — lk V i 2 — l,l >2 tty 2 /r V i 2i y 2 — ] J 


+«■ £ 


n + hot _ „ 

n V 2 - v K V2 _ VV2 n 


»2 




hot I 

V2,V2~V I 


(4) 


where a is the fraction of total 0( 3 P) density, which corre- 
sponds to hot atoms, ky°l v lh and ky, ot r ,_y are the rate coef- 
ficients for excitation and de-excitation of CO 2 molecules, 
respectively, due to collisions with hot atoms, assuming also 
multi-quantum processes. These rate coefficients are not re- 
lated by the detailed balance since hot 0( 3 P) atoms are 
not thermalized. Comparing Eq. (2) which is applied in the 
model used in our study with Eq. (4), one can see that the 
rate coefficient values retrieved in this work and in other at- 
mospheric studies are some sort of effective rate coefficient. 
This may be expressed as 


*£o(z) = C-ife) = d - “(*)) ■ k n-hv 2 

+a(z)-J2 n n k vT-v,v 2 

V 


(5) 


which includes the contribution of hot Of 3 P) atoms. We note 
here that the k\]\"_ v |: , in Eq. (5) relates only to the pump- 
ing term and, therefore, kJ e {L strictly speaking, should not 
be treated as a reaction coefficient of the balanced process 


anymore, which is important for the future studies of this 
problem. In addition to non-thermal oxygen, one can also 
consider other sources of collisional excitation as, for ex- 
ample, collisions with thermal and non-thermal hydrogen, 
electronically-excited atomic oxygen 0( 1 S), charged compo- 
nents, or multi-quantum excitation of CO 2 by thermal oxy- 
gen (Ogibalov, 2000). In any case, the effective collisional 
rate coefficient and the VMR of this “unknown” component 
with respect to Of") may be represented in a way similar to 
Eq. (5). However, in the case of multi-quantum excitation 
y 2 must be replaced with the multi-quantum excitation 
rate, while a(z ) in this approach should be constant with alti- 
tude, which is in contradiction with variable kyj(z) observed 
in the present study. 

Simple calculations show that assuming k Vl - i, v , to be 
fixed at (5.0) x 1 0“ 1 3 [cm 3 s -1 ], the rate coefficients for the 
collisions with hot oxygen, k\° i _ v V2 , should be 10 times 
larger than k m - i iV2 and a(z) should vary linearly from 0% at 
85 km to 2 % at 105 km in order to account for the observed 
kyi{z) increase with altitude. As Kharchenko et al. (2005) 
show, the hot oxygen fraction at these heights is ~ 5 % dur- 
ing daytime. Even though most of the overlapping profiles 
used in this work fall in 18-6 h local time interval, the MLT 
area is illuminated by the sun in at least half of the se- 
lected events (the shadow level at 100 km corresponds to so- 
lar zenith angle Q z — 101°, and the average 6- is 105 + 30°). 
The magnitudes of k„° l v Vl are larger than k V2 - i,y 2 (Shizgal, 
1979) though obtaining the exact values requires complicated 
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quantum-mechanical calculations (V. Kharchenko, personal 
communication, 2011). It is also important to calculate the 
thermalization timescales for hot oxygen atoms. To summa- 
rize: the concept of hot oxygen as an additional pumping 
source for the CC> 2 (i' 2 ) levels is a possible explanation of 
the kyj(z) variation and the discrepancy between laboratory 
measurements and atmospheric estimates of this rate coeffi- 
cient. Of course, the simplified analysis given above is not 
intended to be the ultimate explanation of our result. How- 
ever, it may indicate the direction toward its interpretation. 

5 Summary 

In this work, we discussed a refined methodology for retriev- 
ing ky t- the quenching rate for the V 2 vibrations of CO 2 
molecules by collisions with O atoms. The suggested ap- 
proach differs from those used in the previous studies (see 
the “atmospheric retrieval” values in Table 1), where the kyj 
was retrieved by fitting the measurements of the CO 2 15 pm 
emissions. The method (a) combines space-born radiance 
measurements with the ground-based (i.e. lidar) temperature 
observations, (b) allows retrieving A?vt independently at vari- 
ous altitudes, (c) allows retrieving the parameter y = kyj x O 
(where O is the VMR of atomic oxygen), which is the con- 
tribution of atomic oxygen quenching to the total quenching 
rate of the CO 2 U 2 ) vibrations. Using the y variable instead 
of the kyj in minimizing the difference between measured 
and simulated limb radiance significantly improves the sta- 
bility and accuracy of retrieving. 

The suggested technique is applied to the overlapping 
SABER/TIMED (15 pm radiance, C0 2 and 0( 3 P) VMRs) 
and Fort Collins lidar temperature measurements. The ob- 
tained y values are 1.1 xlO -14 ; 5.7 x 10 -14 ; 2.0 xlO -13 ; 
5.1 x 10 -13 ; 6.4 x 10 -13 [cm 3 s _1 ] at 85; 90; 95; 100; 
105 km, respectively. The average value of the retrieved 
kyj in 85-105 km altitude range is 6.5 + 1.5 x 10“ 12 cm 3 s -1 
that is in excellent agreement with other values esti- 
mated from space-borne observations. However, we also 
observed an altitude dependence of the retrieved kyj, 
which varies from 5.5 + 1.1 x 10 -12 cm 3 s~ [ at 85km to 
7.9 + 1.2 x 10 -12 cm 3 s _1 at 105 km. We link the observed 
variation as well as the discrepancy between the “atmo- 
spheric” and laboratory measurements to a simplification 
in the traditional consideration of C02(t’2)+ 0( 3 P) interac- 
tions in the current non-LTE models, including the model 
applied in this study, where the rate coefficients of the vi- 
brational quenching and excitation are linked by the detailed 
balance relation. We discuss that both problems may be ex- 
plained by interaction of CO 2 molecules with the hot oxy- 
gen atoms, which is missing in the current models. We sub- 
divide the oxygen atoms in the MLT to “normal” and “ex- 
cited” or “hot” groups, with different ky t rate coefficients 
for each of the groups and keep in mind that rate coefficients 
for the former group are balanced, whereas those for the lat- 


ter group do not obey the detailed balance relation. Further 
studies to explain the current results are needed including (a) 
quantum-mechanical calculations for the collisions of CO 2 
with “hot” O atoms and, possibly, other atmospheric compo- 
nents needed to improve the non-LTE model of the IR emis- 
sion formation in the MLT, (b) new comparisons of space 
radiance observations with the lidar temperature measure- 
ments at various locations (including those for polar regions). 
The radiative cooling rate calculations for general circulation 
models and calculations of the emerging IR radiation will ob- 
viously depend on the mechanism of interaction between the 
radiation field and atmosphere, which will be revealed in the 
course of these additional studies. However, one result fol- 
lows already from the present work: for the practical temper- 
ature retrievals from the 15 pm CO 2 atmospheric radiation 
observations, the A'J e g(z) values obtained in this study may 
be recommended. As we show in our discussion, at the cur- 
rent stage of the non-LTE model development these values 
provide better presentation of the CO 2 -O interactions, which 
relies at best fitting of coincidental atmospheric radiance and 
ground based temperature observations. 

Appendix A 


Abbreviations 


ACE FTS 

Fourier-Transform Spectrometer in 
Atmospheric Chemistry Experiment 

ALI-ARMS 

Accelerated Lambda Iterations for Atmo- 
spheric Radiation and Molecular Spect 

ASTRO-SPAS 

Astronomical Shuttle-Pallet Satellite 

CRISTA 

Cryogenic Infrared Spectrometers and 
Telescopes for the Atmosphere, on board 
the ASTRO-SPAS satellite 

EUV 

Extreme Ultraviolet 

LTE 

Local Thermodynamic Equilibrium 

MGS TES 

Thermal Emission Spectrometer onboard 
Mars Global Surveyor satellite 

MIPAS 

Michelson Interferometer for Passive At- 
mospheric Sounding 

MLT 

Mesosphere/Lower Thermosphere 

SABER 

Sounding of the Atmosphere using 
Broadband Emission Radiometry 

TIMED 

Thermosphere Ionosphere Mesosphere 
Energetics and Dynamics 

VMR 

Volume Mixing Ratio 

WACCM 

Whole Atmosphere Community Climate 
Model 
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